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Abstract 

The common rules of classical nonequilibrium thermodynamics do not allow 
an Onsager coefficient like the viscosity to depend on the shear rate. Such a 
dependence is experimentally well documented, however. In this paper it is 
shown, using nonequilibrium thermodynamics alone, how it is possible that 
the effective viscosity coefficient of an isotropic colloidal system depends on 
both the shear rate and the frequency. 

I. INTRODUCTION 

Much is known about the hydrodynamic description of colloidal systems. This is de- 
scribed in a number of well written books on the subject [1,2,3]. In these books a multitude 
of details about the interesting behavior of suspended asymmetric colloidal particles in 
shear flow is treated. What is not being done is to treat these systems in the context of 
nonequilibrium thermodynamics. That such a treatment is possible becomes clear reading 
the monograph by De Gennes [4,5] on The Physics of Liquid Crystals. Rather than going 
into detail regarding the motion of individual particles he discusses, both the equilibrium 
and the nonequilibrium properties of these materials, using thermodynamics. In this paper 
we will apply nonequilibrium thermodynamics to these colloidal systems. As it is not our 
intention to duplicate all the well described knowledge about the collective behavior of these 
systems, we will focus on one aspect which has remained puzzling from the viewpoint of 
nonequilibrium thermodynamics. 

It is well-known that the viscosity in these systems depends on the shear rate. In view 
of the fact that the shear rate is not a thermodynamic state variable, like, for instance, the 
temperature, classical nonequilibrium thermodynamics, as explained in the classic mono- 
graph by de Groot and Mazur [6], does not allow such a dependence. Books on colloidal 
hydrodynamics explain this dependence in terms of the behavior of asymmetric suspended 
particles in shear flow. The question we pose ourselves in this paper is: explain this de- 
pendence using nonequilibrium thermodynamics alone. Interestingly enough Prigogine and 
Mazur [7] already mentioned this problem as one of the three possible examples of their 
method of internal variables. The orientation of the suspended particles was suggested as 
such an internal variable. In de Groot and Mazur's book [6] the other two examples are 
treated in detail but this one was not pursued. It had to wait for a better understanding 
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of these systems. In this paper we will show that it is now possible to answer this question 
for a suspension of asymmetric colloidal particles. In addition to answering a question of 
principle, the treatment is simpler and more direct than the current explanations. It is also 
not restricted to dilute suspensions. 

In the second section we give the Gibbs relation for an isotropic system with asymmetric 
colloidal particles. The particles are assumed to have a symmetry axis. A new term appears 
which gives the entropy change due to a change in a tensorial order parameter Q. This 
order parameter characterizes the alignment of the particles [4,5]. This is precisely the 
new internal variable Prigogine and Mazur [7] needed. In the third section the resulting 
entropy production is presented. This leads to linear laws in which the viscous pressure 
tensor and the time rate of change of the order parameter are expressed in their conjugate 
thermodynamic forces. In the fourth section these laws are solved for the special case of 
stationary shear flow. An expression is found for the viscous pressure in terms of the shear 
gradient. Their ratio defines an effective viscosity. This effective viscosity is found to be a 
function of the shear. Crucial to obtain this result is the fact that the Onsager coefficients 
in general depend on the order parameter Q which is a state variable. This answers the 
principle question. We proceed in chapter five with a analysis of time dependent shear flow. 
It is found that this results in a frequency dependent effective shear viscosity. In the last 
section we derive a simple expression for the dependence of one of the Onsager coefficients 
on the order parameter to make this point clear beyond doubt. 



In an isotropic suspension of particles, which are cylindrically symmetric, a tensorial or- 
der parameter Q appears as a new thermodynamic variables [4,5]. Q is a traceless symmetric 
tensor, which is defined as the ensemble average of (aa — 1/3) , where a is a unit vector along 
the symmetry axis of one of the particles and 1 is the unit tensor. In an isotropic phase Q 
is, in the absence of any directing field, zero. In, for instance, a magnetic field Q is nonzero. 
Changes of Q lead to changes of the entropy. It is therefore important to take the order 
parameter along in the description of the system. The total differential of the entropy is 
given by the following Gibbs relation 



Here T is the temperature, s v = ps the entropy density, u v = pu the internal energy 
density, p the chemical potential of the solvent, p the solvent particle number density, p s 
the chemical potential of the suspended particles and p s the suspended particle number 
density. The densities with a subscript v, in addition to p and p s , are per unit of volume, 
w is the variable conjugate to Q and is given by 



II. THE GIBBS RELATION 



Tds v = du v — pdp — p s dp s — w : o?Q 



(1) 




(2) 



Alternatively one may use the relation to the free energy /, 



V 




(3) 
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Furthermore : between two matrices signifies a double contraction. De Gennes [4 (page 
48), 5] writes as explicit expression for the free energy density 

/„ (T, p, p„ Q) = f v (T, p, p s ) + l -A (T, p, p s ) tr (Q.Q) +±B (T, p, p s ) ir (Q.Q.Q) + O (Q 4 ) 

(4) 

where O is the order symbol. Here tr signifies the trace of a matrix. Note that tr (Q.Q) = 
Q : Q. The resulting traceless symmetric conjugate variable is 

w = A (T, p, p s ) Q+B (T, p, p s ) CpJ + O (Q 3 ) (5) 

o 

where Q.Q is the traceless symmetric part of Q.Q. This relation may be inverted to give 
the order parameter in terms of its conjugate variable 

Q = A' 1 (T, p, p s ) w-A- 3 (T, p, p s ) B (T, p, p s ) w^w + O (w 3 ) (6) 
The typical size of A is of the order of k B Tp s where ks is Boltzmann's constant. 



III. THE ENTROPY PRODUCTION 

Without the Q term the Gibbs relation is equivalent to the one used in de Groot and 
Mazur [6]. Note that they use densities per unit of mass while we use densities per unit of 
volume like de Gennes [4,5] is doing. The calculation of the entropy production, using the 
balance equations, is essentially identical to the one they give. The only difference is the 
appearance of a term due to the order parameter Q. In order to describe the contributions 
of the order parameter to the dynamical behavior of the system we may restrict ourselves, 
using the Curie principle, to contributions of force-flux pairs which are symmetric traceless 
tensors. The rate of entropy production density, as a function of position r and time t, due 
to these contributions becomes 

dQ (r,t) 



Ta tens (r, t) = - II (r , t) :gradv (r, t) -w (r, t) : ' (7) 



here II is the symmetric traceless part of the viscous pressure tensor and gradv the sym- 
metric traceless part of the velocity gradient. 
The resulting linear laws are 



II (r, t) = —2r] gradv (r, t) —£w (r, t) 
dQ(r,t) 



dt 



I gradv (r, t) -£ Q w(r,t) (8) 



Here we used that, gradv (r, t) is odd and w (r, t) even for time reversal, so that the matrix 
of Onsager coefficients is antisymmetric. The Onsager coefficients are in general functions 
of the thermodynamic state variables, so that we have 
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rj{T,p,p s ,Q) , £(T,p,p s ,Q) and i Q (T,p,p s ,Q) 



(9) 



In line with the assumption of linearity, they are not allowed to be functions of gradv (r, t) 
and dQ (r, t) jdt. The diagonal coefficient r\ is the viscosity and the diagonal coefficient £q is 
a typical rotational mobility of the suspended particles divided by p s . The cross coefficient 
I is dimensionless and is expected to be a positive constant of the order unity. 



IV. STATIONARY ELONGATION AL FLOW 



The above equations may easily be solved for stationary elongational flow 



v(r,t) = - 7 (y,x, 0) 



The resulting velocity gradient is 



gradv (r, t)= gradv (r, t) — - 7 



(° 1 °\ 
1 

\o o o; 



Using the second linear law, eq.([8]b), results in 



w (r, t) = — gradv (r, t)- 



2L 



1 
- | 1 




Substitution in eq.(BI) gives 



Q 



2A£, 



7 



/0 1 0\ 

1 + O (7 2 ) 
\0 / 



n( r ,t)= 



rj + 



lie 



/0 1 
1 
\0 



(10) 



(11) 



(12) 



(13) 



(14) 



The first linear law, eq.(|8]a), then yields 

/ „■>. x / 1 

7 10 0]= -Tj eff 7 
\0 

Writing the dependence on the thermodynamic state variables explicitly, the effective vis- 
cosity becomes 

1,2 (T,p,p s ,Q) 



Veff (T, p, ps, Q) = V ( T i Pi Ps, Q) + 



2£ Q (T,p,p„Q) 



(15) 



In view of the above relation of the order parameter with the rate of elongation, one may 
write this relation as 



Veff (T, p, p s , 7) = ?7 ( T , P, Ps, + 



(t, p, p s 



7 



2£ Q (T,p,p s ,7 

The effective viscosity depends therefore on the rate of elongation. 



(16) 
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V. OSCILLATORY ELONGATIONAL FLOW 



The above equations may also easily be solved in the case of oscillatory elongational flow. 
In that case the velocity field, its symmetric traceless gradient, the order parameter and its 
conjugate variable all become proportional to exp(— iut). We now neglect the nonlinear 
terms in eq.(§) and the Q dependence of the transport coefficients in order to avoid the 
appearance of higher order harmonics. Though this is an interesting phenomenon, it is 
not our present concern. In the second linear law we must now replace dQ (r, t) / dt by 
— iuQ (r, t). In the solution this implies that 

t - 

Q =A-'w = — gradv (17) 

A£ Q - iuj 

The resulting complex frequency dependent effective viscosity then becomes 

Veff(T,p,p s ,uj)=r](T,p,p s )+ Z (Tip,p 8 ) 

(T, p, Ps) (1 - iujt (T, p, p s )) 

with a relaxation time 

r (T, p, p s ) = — — — — (19) 

A(T,p,p s )£ Q (T,p,p 8 ) 

This relaxation time is the reorientation time of the suspended particles due to rotational 
diffusion, and is therefore of the order of the rotational diffusion coefficient. It follows that 
lot is the relevant Peclet number for rotational diffusion in an oscillating velocity field. For 
low frequencies the effective viscosity is enhanced by the coupling to the order parameter 

Veff (T, p, p., u = 0) = rj (T, p, p.) + f p { f T P,Ps \ (20) 

M Q {i-iPi Ps) 



This is the same result as found above, eq. (|15|) , for the stationary case if one neglects the 
dependence of the coefficient on the order parameter. For high frequencies one finds 

Veff(T,p,p s ,uj = oo) = T](T,p,p s ) (21) 

Writing the complex viscosity as the sum of a real and imaginary part, r) e ff = f]' e ff + Z7 7e//' 
one has 

Veff( T iPiPs,u)=ri(T,p,p s )+ Mt~T^2 (rr ^ ( 22 ) 

" 2£ Q (T,p,p a ) (1 + uj 2 t 2 (T,p,p s )) 

and 

// frp \ u£ 2 (T,p,p s ) 

Veff {T,PiPs,U) = 23 

n 2£ Q (T,p,p s )(l + ^ 2 r 2 (T,p,p s )) 

This last expression gives the elastic contribution to the complex viscosity. 
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VI. THE ONSAGER COEFFICIENTS 



Nonequilibrium thermodynamics does not derive expressions for the Onsager coefficients. 
Like the thermodynamic derivatives of the entropy as for instance A (T,p,p s ), they have to 
be either measured or calculated from a description in terms of the motion of the separate 
particles. Above we have given estimates of the typical order of magnitude of the vari- 
ous coefficients on the basis of the physical meaning of the coefficients and dimensional 
arguments. 

As the shear dependence of the effective viscosity is the property we set out to find, we 
will pursue this element a bit further. One may write 

V = 77o (1 + [77] 0) (24) 

where = Vp s , with V the volume of the monodisperse suspended particles, [77] is the 
so-called intrinsic viscosity. Each suspended particle contributes to the viscous dissipation. 
For low densities the intrinsic viscosity is independent of their concentration. Einstein [8] 
found that the intrinsic viscosity for spheres was equal to 2.5. For rods the contribution 
depends on the orientation a of the rod. In order to get the intrinsic viscosity one must 
then first calculate the contribution of a rod to the viscous pressure tensor as a function of 
its orientation and then average over the orientations [2]. 

The distribution over orientations may be calculated using the rotational diffusion equa- 
tion for this distribution. For the stationary state distribution corresponding to a given 
order parameter Q one has 

'< a > = 1 (— I) : Q+ i (25) 

For spheroids, which are useful as a model system, the solution of the velocity field has 
been given by Jeffery [9] . In the monograph by van de Ven [2] the resulting contribution to 
the pressure tensor is given. In the appendix we average those expressions over the above 
orientation distribution. This results in 



u 7 



4 (2r e 2 - (1 - 2r 2 ) A e ) A£ Q 

p, 

(26) 



26r 2 e - 24r 2 A e - 15A e 6 



(2r 2 - 3A e ) (2r 2 - (2r 2 + 1) A e ) (r 2 + 1) (3A e - 2 

Here the aspect ratio r e = a/b is the ratio of the diameter a along the symmetry axis and 
the diameter b normal to the symmetry axis. Furthermore 



r 2 r e arccosh r e 

A e = — tjtt for r e > 1 

r-2-1 (r 2 -l) 3/2 



e 



r e arccosr e r 2 „ . . 

A e = -| e — for r e < 1 27 

(1-rlf 2 l-r e 2 K J 

Note that 7 /AIq is the rotational Peclet number so that eq.(|26|) is an expansion in the 
Peclet number. In the limit of a sphere, r e — > 1, one may verify that the expression reduces 
to 2.5. In that case the shear dependence of the effective viscosity disappears. Prolate 
spheroids, r e > 1, are shear thickening and oblate spheroids, r e < 1, are shear thinning for 
small Peclet numbers in elongational flow [10]. 
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VII. DISCUSSION AND CONCLUSIONS 



The main objective was to show that a straightforward application of nonequilibrium 
thermodynamics leads to a shear dependent viscosity. We were able to derive a general 
nonlinear expression for this dependence. The important reason for this is the coupling 
to the orientation of the suspended colloidal particles. This orientation is described by a 
tensorial order parameter Q as thermodynamic state variable. The fact that the Onsager 
coefficients depend on this state variable then results in a shear dependent viscosity. Of 
course the critical reader could have doubted such a dependance of the Onsager coefficients 
on Q. In order to defend us against this critique, we derive an explicit formula for small 
Peclet numbers, 7 /AIq. This gives a linear relation for this dependence. For larger Peclet 
numbers the nonlinear relation between Q and w must be taken into account. This will lead 
to more realistic and interesting predictions. One interesting aspect is that the powers of Q 
lead to shear field with a different "direction". As the principle question, why the viscosity 
depends on the shear, has been answered, we will not here pursue these points. 

The frequency dependence is due to a single exponential decay of the order parameter. 
For low frequencies the viscosity is enhanced by the coupling to the orientation of the colloidal 
particles. For higher frequencies the orientation lags behind. This leads to thinning of the 
fluid and a reversible elastic response. One can extend the analysis, within the general 
frame of nonequilibrium thermodynamics, by introducing the whole angular distribution of 
the particles as an internal thermodynamic variable [6]. The details for such an analysis 
have recently been worked out by Mazur [11]. This results in a whole spectrum of relaxation 
times and consequently gives a more complex frequency dependence in the transition from 
low frequency to high frequency behavior. 

In our analysis we do not need to analyse the orientational distribution function in 
any detail. The knowledge of the order parameter is enough. This makes our analysis much 
simpler that the usual procedure, where the rotational diffusion equation for this distribution 
function is solved numerically. As we said above [11], the orientational distribution function 
can be introduced as an internal thermodynamic variable. In addition to a richer frequency 
dependence this will also lead to a more complex shear dependence. 

The low frequency and the low shear behavior of the viscosity has a very similar non- 
analytic dependance on these parameters [12]. This is a general phenomenon which is not 
specifically due to the asymmetry, or for that matter the presence, of the suspended particles. 
It falls outside the scope of the present article. 

In the field of socalled extended irreversible thermodynamics [13], one assumes that 
fluxes, fluxes of fluxes, etc., are also thermodynamic state variables. This assumption dis- 
tinguishes it from classical irreversible thermodynamics, where this is not considered to be 
a proper assumption. The viscous pressure is such a flux and is, according to extended 
irreversible thermodynamics, a state variable. On the basis of this assumption the viscosity 
can be a function of the viscous pressure. As a result, the viscosity then becomes a function 
of the shear. The assumption, so to say, defines the problem away. We will not try to inval- 
idate their assumption. We only want to stress that the present paper shows that such an 
assumption is not needed to obtain a shear dependent viscosity. As such their assumption 
seems to be one to many. 
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VIII. APPENDIX 

We use the results of Jeffery [9] as presented in van de Ven's monograph [2], see in 
particular pages 229-231. The shear field he uses is V\ = — and = Gx<i- In 
our notation x — x^, y — x%, z — x%, (j) — c and 7= G. The symmetric traceless part of 
this simple shear field is identical to the one we used above. In view of the linearity of the 
problem we may therefore use the symmetric traceless part of the contribution to the viscous 
pressure tensor van de Ven uses: 



n = n r 



57TV70 

IT 1 



<M> 



(28) 



where V = 4irb 3 r e /3 is the volume of the particle. IIq = —2t]q gradv (r,t) is the viscous 
pressure one would have in the absense of the suspended particles. He gives the tensor M 
as M' in eq. (3.280) in a coordinate frame attached to the particle. Using the distribution 
given above, eq.(|25"D, we have 



(M) 



15 
8^ 



M (an aa 



dQ 



: Q + 



47T 



M fa) dQ 



(29) 



where dQ = sm9d9d(fi with < 9 < n and < < 2n. The integral between square 
brackets gives a constant times the four index symmetric traceless unit tensor A. The 
matrix elements of this tensor are given by 



A 



ijkl 



(30) 



The constant can be calculated by taking i 
results in 



I, j = k and summing over i and j. This 



(M) 



16n 



M (a) : aa 



dQ 



Q 



l 

47T 



M fa) dQ 



(31) 



The contraction of two matrices is independent of the reference frame and we may therefore 
use the reference frame attached to the particle. In this reference frame (aa— 1/3) = 
(2/3) eiei — (1/3) e2e2 — (1/3) 6363 in terms of the unit vectors in this frame. This gives 



(M) 



16n 



(2M' U -M^ 2 (8,<t,)-M! a (8,<t,))dn 



M 



dQ 



(32) 



Using eq. (3.280) we see that M' 22 and M33 integrate to zero while M' 1X is independent of 
direction. In the second integral we can use eqs.(3.282)-(3.284) to obtain the only nonzero 
contribution. The result is 



S 



1 b 3 r / 4 1 2 

<m) = -m; iQ + ^ (- Ci + -c 2 + -c 3] - 



/O 1 
1 
\0 



Substitution of Q then yields 
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/0 1 
1 

Vo o o 



Substitution into eq.(|28"D and comparing with eq. (p4|) gives for the intrinsic viscosity 



3M' 



Ai 2b 3 r e 15 



11 + — Ci + + =-C 3 



1 

—i 

3 



2 

-( 

3 



Substitution of eqs.(3.280) and (3.284) finally gives eq.©. 
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